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Abstract. Magnetic fields can be created in stably stratified (non-convective) layers in a differentially rotating 
star. A magnetic instability in the toroidal field (wound up by differential rotation) replaces the role of convection 
in closing the field amplification loop. Tayler instability is likely to be the most relevant magnetic instability. A 
dynamo model is developed from these ingredients, and applied to the problem of angular momentum transport 
in stellar interiors. It produces a predominantly horizontal field. This dynamo process is found to be more effective 
in transporting angular momentum than the known hydrodynamic mechanisms. It might account for the observed 
pattern of rotation in the solar core. 
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1. Introduction 

Magnetic field generation by differential rotation is usu- 
ally regarded as a process operating in the convective 
zones of stars. A 'convective stellar dynamo' (e.g. Parker 
1979) is then a field amplification process in which dif- 
ferential rotation stretches field lines into a toroidal field 
(i.e. an azimuthal field, with respect to the rotation axis). 
Convective fluid displacements create a new poloidal field 
component by putting kinks into the toroidal field lines. 
The stretching of these kinks produces a new toroidal 
component, amplifying the existing toroidal field and thus 
'closing the dynamo loop'. The process is widely believed 
to be responsible for the magnetic fields seen at the sur- 
face of stars with convective envelopes like the Sun, and 
the study of this kind of convective field generation has 
created a literature of vast proportions over the past 5 
decades. 

Convection, or other imposed velocity fields like waves 
or shear turbulence, are not really necessary, however, for 
a dynamo process to operate. An example of a magnetic 
field produced by differential rotation without the assis- 
tance of an imposed small scale velocity field is the field 
produced in accretion disks (Hawley et al. 1996). In this 
kind of dynamo process, the differential rotation produces 
a small scale magnetic field on its own: the role of an 
imposed velocity field in generating a new the poloidal 
field component is replaced by an instability in the mag- 
netic field itself. The magnetic instability operates on the 
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toroidal field that is produced by winding up of the radial 
field component, distorts it, and creates new radial field 
components. The instability can either be of the Velikhov- 
Chandrasekhar-Balbus-Hawley type (hereafter BH insta- 
bility, see Velikhov 1959; Chandrasekhar 1960; Balbus & 
Hawley 1991), or a buoyancy- driven instability (Parker 
1966), though the properties of the resulting field may de- 
pend on which of these instabilities is the most important. 

In the same way, the generation of a magnetic field in 
a star requires only one essential ingredient: a sufficiently 
powerful differential rotation. The recreation of poloidal 
field components which is needed to close the dynamo loop 
can be achieved by an instability in the toroidal field, of 
which there exists a variety including those (BH and buoy- 
ancy instabilities) that operate in disks. Such instabilities 
do not require the presence of convection, and can take 
place also in the stably stratified interior of a star. 

The instabilities of a predominantly toroidal field in 
stably stratified regions in a star have been summarized 
from the existing literature in a previous paper (Spruit 
1999). I concluded there that the first instability likely to 
set in is a pinch-type instability. The essential properties 
of this instability in the stellar context were established 
by R.J. Tayler (Tayler 1973; Markey & Tayler 1973, 1974; 
see also Tayler 1957; Goossens et al. 1981). What makes it 
of particular importance are the absence of a threshold for 
instability (at least in the absence of viscous damping and 
magnetic diffusion; more about this below), and the short 
growth time, of the order of the Alfven crossing time. It 
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can can operate under conditions where BH and buoyancy 
instabilities are suppressed by the stable stratification. 

As an aside I note here that the solar cycle, generally 
considered as the classical case of a convectivc dynamo 
process, is probably not driven by convective turbulence 
at all. Properties of active regions such as the heliographic 
latitude and the time scale at which they emerge, and the 
value of the tilt of their axes, point to field eruption from 
an ordered toroidal field at the base of the convection zone. 
The energy density in this toroidal field inferred from the 
observations is about 100 times stronger than can be ex- 
plained by convective turbulence (D'Silva & Choudhuri 
93; D'Silva and Howard 1993). The relevant processes in 
such a strong field are the winding-up by differential rota- 
tion, and buoyant instability (Parker 1966) of the field it- 
self (Caligari et al. 1998 and references therein). This view 
is much closer to semiempirical models of the cycle like 
those of Leighton (1969), the development of which was 
eclipsed by the mathematically more interesting turbulent 
dynamo models. Leighton's model is also much closer to 
the observations, since it includes the role of active region 
eruption as an essential ingredient in the operation of the 
dynamo process. This contrasts with the convective dy- 
namo view, in which active regions are only a secondary 
manifestation of turbulence. 

1.1. Energetics of magnetic field generation 

The difi^erential rotation of a convectively stable stellar 
interior is 'limited resource'. Whereas a dynamo operat- 
ing on the differential rotation in a convective zone feeds, 
indirectly, on the continuous energy flux of the star, a dy- 
namo process in stable layers derives its energy from the 
rotation of the star. This is a very small amount of energy 
compared with the luminous energy, integrated over the 
star's life time. The energy budget of a dynamo process 
in a stable zone is thus limited. 

The rotational energy is available for a dynamo process 
only though the variation of the rotation rate through the 
star. Such a gradient can result from a spin-down torque 
on the star caused by angular momentum loss by a stellar 
wind, or from changes in the distribution of the moment 
of inertia of the star as a result of its evolution. 

1.2. Angular momentum transport 

The main interest of a magnetic field generated in the 
stable layers lies in its ability to exert internal torques. If 
such torques can operate on time scales short compared 
with the spindown time scale or the stellar evolution time 
scale, respectively (depending which of these processes is 
the main source of differential rotation) , they can strongly 
affect the degree of differential rotation that can survive 
in a star. In some cases, they may maintain a state of 
nearly uniform rotation. For example, the usual interpre- 
tation of the rotation of white dwarfs and neutron stars 
as a remnant of the initial angular momentum of their 



main-sequence progenitors would require a very large de- 
gree of differential rotation between core and envelope of 
these stars during their evolution on the giant and AGB 
branches. Spruit & Phinney (1998) have argued that, in 
the absence of more detailed knowledge, angular momen- 
tum processes are more likely to be either too ineffective 
or too effective, compared with the intermediate effective- 
ness needed to explain the rotation rates of white dwarfs 
and neutron stars. If the process is ineffective, leading 
to evolution at approximate conservation of angular mo- 
mentum, the cores of stars on the AGB are very rapidly 
spinning, and much too rapidly rotating white dwarfs and 
neutron stars are the result. On the other hand, effective 
coupling leads to nearly uniform rotation and too slow 
rotation of the remnants. In the latter case, other mecha- 
nisms must contribute to the rotation of white dwarfs and 
neutron stars, such as birth kicks for neutron stars (Spruit 
& Phinney 1998), and nonaxisymmetric mass loss on the 
AGB for white dwarfs (Spruit 1998). 

All known non-magnetic angular momentum transport 
processes in stars are weak except at rotation near the 
breakup rate. If only these known processes are included, 
pre-supernova cores rotate much too fast to explain the 
majority of supernovae (Heger et al. 2000), though they 
may rotate at the right rate for the coUapsar model of 
gamma-ray bursts (McFadyen, Woosley & Heger 2001 and 
references therein). 

Microscopic viscosity is negligible in stars as a source 
of internal torque. Hydrodynamic processes such as circu- 
lation (Zahn 1992; Maeder and Zahn 1998) and hydrody- 
namic instabilities arc potentially important. In practice, 
however, the torrjues transmitted by most of these pro- 
cesses scale as (Q/N)^ (where fl is the rotation rate and 
N the buoyancy frequency) and hence are rather ineffec- 
tive except in stars rotating close to their breakup speed, 
and can not explain the low rotation rate of the solar core 
(Spruit et al. 1983). An exception is angular momentum 
transport by internal waves generated by pressure fluctu- 
ations in a nearby convective envelope (Zahn et al. 1997), 
which scales as a lower powcir of fl. Even these waves, how- 
ever, appear insufflcient to explain the near- uniform rota- 
tion observed in the solar interior (Talon & Zahn 1998; 
Kumar et al. 1999). 

1.3. Evidence from the internal rotation of the Sun 

The very small degree of differential rotation in the core 
of the Sun, and the small difference in rotation rate be- 
tween the core and the convective envelope (Schou et al. 
1998; Charbonneau et al. 1999) requires the presence a 
process with two special properties. First, the envelope is 
continuously spun down by the solar wind torque, whose 
magnitude can be determined from in-situ measurements 
in the wind (Pizzo et al. 1983). The absence of an increase 
of rotation with depth in the core, and the near-equality 
of the rotation in the core to the average rotation in the 
convective envelope, require an efflcient angular momen- 
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turn transport process in the core. Secondly, there is no 
noticeable difference of the core's rotation with latitude, in 
contrast with the convective envelope, which rotates some 
40% faster at the equator than at the poles. An isotropic 
angular momentum transport process coupling the core 
to the differentially rotating envelope would lead to rota- 
tion with roughly equally strong vertical and horizontal 
gradients. This is contrary to observation: at the transi- 
tion between core and envelope (the so-called tachocline, 
Spiegel & Zahn 1992) the radial gradient is much stronger 
than the horizontal gradient. A prospective angular mo- 
mentum transport process explaining this must bo even 
more effective at reducing the gradient of rotation on hor- 
izontal surfaces than it is in reducing the radial gradient. 

1.4. Magnetic angular momentum transport 

Magnetic fields are quite effective at transporting angu- 
lar momentum, even at strengths much below the values 
observed at the surfaces of magnetically active stars or 
the magnetic A stars (Mestel 1953). The Maxwell stresses 
transport angular momentum in the radial direction at 
a rate proportional to B^B^, and in latitude at a rate 
proportional to BgB^. If the observed internal rotation 
pattern in the solar core is due to magnetic stresses, the 
weak gradient in latitude compared with the radial gradi- 
ent then indicates that Bg Br in the tachocline. This 
fact, in combination with a sufficient amplitude of the 
magnetic field, are to be explained by a prospective theory 
for magnetic angular momentum transport in the core of 
the Sun. 

In the following I show that a differential rotation- 
driven dynamo operating with magnetic instabilities can 

do just this, and present a quantitative estimate of the 
amplitude, length scales and time scales of the resulting 
magnetic field. 

2. The dynamo process 

I assume that the star's rotation rate is a function of 
the radial coordinate r only ('shellular' rotation, cf. Zahn 
1992). This is done for convenience only. It is quite con- 
ceivable that latitude-dependent differential rotation, for 
example due to a nearby convective zone, would have 
somewhat different effects. 

Another assumption is that the initial magnetic field is 
sufficiently weak, so that magnetic forces can be neglected 
initially. Thus, I ignore the possibility that the star started 
with a strong magnetic field, such as those of the magnetic 
A-stars or the magnetic white dwarfs. 

On account of its weakness, the radial component Br of 
the field is wound up by differential rotation. After only a 
few differential turns, the resulting field is predominantly 
azimuthal, B^ J> Br, and its strength increases linearly 
with time, until it becomes unstable. As long as the mag- 
netic field is weak, the energy to be gained from the field 
by any instability is small, and instability is possible only 
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for displacements that avoid doing work against power- 
ful restoring forces: the gas pressure and the buoyancy 
force. The first instability to set in (for a review of the 
relevant instabilities see Spruit 1999) is a nonaxisymmet- 
ric (typically m = 1) interchange in which the displace- 
ments are nearly incompressive (div^ ^ ^/r, thus avoid- 
ing work against the gas pressure), and nearly horizontal 
(along equipotential surfaces), avoiding work against the 
buoyancy force. 

In the absence of diffusive damping effects such as mag- 
netic diffusion and viscosity, an arbitrarily weak azimuthal 
field is unstable (Tayler 1973; Pitts & Tayler 1986). The 
form of the instability is related to kink instabilities in a 
linear pinch. For a purely azimuthal field, the instability 
is a local one: unstable displacements exist as soon as a 
local instability condition is satisfied at any point in the 
star (Tayler 1973). Under the present astrophysical condi- 
tions of a strongly stabilizing buoyancy force, I will refer 
to pinch-type instability of an azimuthal field under strati- 
fied conditions as they exist in a stellar interior as a Tayler 
instability for short. 

2.1. time scales 

There are a number of different time scales in the problem. 
The fastest time scale is the sound travel time over, say, a 
pressure scale height H, tg = {g/H)~^^'^, where g is the ac- 
celeration of gravity. Next slower is the buoyancy (Brunt- 
Vaisala-) frequency, N = tJ^(Vad — V)^/^, where V and 
Vad are the usual logarithmic temperature gradient and 
its adiabatic value (e.g. Kippenhahn and Weigert 1990). 
This frequency is of the same order as the inverse sound 
time except in and close to a convective zone. Since we are 
dealing here exclusively with stable, radiatively stratified 
zones. A'' > is a large frequency. 

Next slower is the rotation rate 17 of the star. Only 
for rotation rates near breakup can fl approach N. The 
lowest frequency is the Alfven frequency, 

'^^ - (4^p)V2^- 

Thus, I will assume the ordering 

WA < n < A^. (2) 

This limit results in significant simplifications and is the 
most relevant one for the processes considered here. Under 
these conditions, the characteristic growth rate of the in- 
stability is of the order 

a = ujl/n {ujA < ^) (3) 

(Pitts & Tayler 1986; Spruit 1999). For slow rotation {n < 
wa) on the other hand, the typical growth rate is wa- The 
reduction by the factor LUj^/fl is due to the strong Coriolis 
force when rotation is rapid. This is characteristic of all 
interchange-type instabilities in a rotating system. 

Finally there is a time scale associated with the differ- 
ential rotation: 

t^^ = {rdrn)-\ (4) 
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2.2. Unstable displacements and instability condition 

The azimuthal field can be seen as consisting of stacks 
of loops concentric with the rotation axis. The instability 
causes loops to slip sideways with respect to each other by 
motions close to equipotential surfaces. The force driving 
the instability is the magnetic pressure of the azimuthal 
field; the loops in a stack press on each other with a pres- 
sure which can be released by sideways displacements 
in a process similar to the accidental slippage of disks in 
the spinal columns of bipedal vertebrates]^. 

A heuristic argument shows how the strong buoyancy 
force (N 3> l^ja) quantitatively constrains the unstable 
displacements. As in the case of hydrodynamic instabili- 
ties, the effect of the buoyancy force is reduced by ther- 
mal diffusion (Zahn 1974, 1983). On small length scales, 
thermal diffusion reduces the stabilizing temperature per- 
turbations caused by the thermal stratification, lowering 
the effective buoyancy frequency. This effect is important 
if the buoyancy is due mostly to the temperature gradient, 
but not if the buoyancy is due mostly to the gradient in 
composition. Thus we have to distinguish these two cases. 
Consider first the case when the buoyancy is due to a 
composition gradient. The effective buoyancy frequency is 
then the same on all length scales. 

2.3. Instability without thermal diffusion 

In this subsection, I ignore the thermal diffusion, but take 
into account the diffusion of the magnetic field by Ohmic 
resistivity. For a dilute (and non-degenerate) plasma, the 
magnetic diffusivity is given by Spitzer's value 



V 



= 7 - lO^lnAT^^/^. 



(5) 



The growth rate of the instability in the absence of 
constraints like rotation is a ~ loa- Thus, the kinetic en- 
ergy released by an unstable displacement of amplitude f 
is ^Wa'^^ (p^^ u'^it; of mass). This energy is supplied by 
the field configuration. In order to avoid wasting this en- 
ergy by doing work against the stable stratification, the 
unstable displacements must be nearly along equipotential 
surfaces, ^h, where = \{^9,^<p)\ is the horizontal 

displacement. If Ir and /h are the radial and horizontal 
length scales of the displacement, the condition that ^ be 
nearly incompressive implies that Ih/lr ~ Ch/Cr ^ 1- For 
such displacements, the work done against the buoyancy 
force is (per unit mass) ^£,'^N'^{lr/lh)'^- Hence the unstable 
displacements have to satisfy Ir/lh < uja/N. 



^ This heuristic picture of the instability applies close to the 
axis, it suggests that unstable displacements cross the axis. The 
instability is actually a local one. If the instability condition 
is satisfied at any point (r, 9) in a meridional plane, there is 
an unstable displacement confined to an azimuthal annulus 
of some width around this point (appendix in Tayler 1973). 
The region near the axis is always unstable; at other latitudes 
the stability depends on the strength of the gradient of field 
strength with latitude. 



In a star, the largest length scale available is of order 
r, so we must have^ Ih < r. So the radial length scale has 
to satisfy 



Ir < toja/N. 



(6) 



On the other hand, if the radial length scale is too small, 
field perturbations cannot grow because the magnetic dif- 
fusivity smoothes them out too fast. This limit is found 
by equating the magnetic diffusion time scale on the 
length scale 1^ to the growth time scale of the insta- 
bility. In presence of a strong Coriohs force, uja, the 
intrinsic growth rate of the instability (in the absence of 
stratification and magnetic diffusion) is cr Ld\/D,. Hence 
the radial length scale has to satisfy 



ll > rin/Lu\. 

Combining the two limits yields 



n 



(7) 



(8) 



This is, for the case when thermal diffusion can be ne- 
glected, and up to a numerical factor of order unity, the 
correct instability condition as derived from Acheson's 
(1978) dispersion relation for azimuthal magnetic fields 
(Spruit 1999). 

2.4. Effect of thermal diffusion 

Consider next the case when composition gradients can 
be neglected, so the stabilizing stratification is due en- 
tirely to the entropy gradient. Since heat is transported 
by photons, while viscosity and magnetic fields diffuse by 
Coulomb interactions, the Prandtl number Pr= v/k and 
the diffusivity ratio s = tj/k are very small in a stellar 
interior (of the order 10~^-10~^). Here v is the viscos- 
ity, K = 16aT^ / {3kb.P^Cp) is the thermal diffusivity, and rj 
the magnetic diffusivity. The viscosity is of the same order 
but somewhat smaller than rj so that the main mechanism 
damping the instability is magnetic diffusion. I ignore vis- 
cous damping in the following. 

If unstable displacements take place on a radial length 
scale / and time scale r, their temperature fluctuations 
diffuse away on the time scale tt — l"^ / k. On the time scale 
T they are therefore reduced by a factor f = {t/tt + 1). 
This can be taken into account by introducing an effective 
thermal buoyancy frequency N^. 



(9) 

An exception would be Z = 1 displacements with ~ 0, 
corresponding to pure rotations of adjacent concentric spheres 
with respect to each other. The radial length scale could be 
arbitrary for such displacements without doing work against 
buoyancy. Since the field configuration is never unstable at all 
latitudes, however, such displacements probably are stable be- 
cause the energy gained from the magnetic field at one latitude 
is lost by work against the field at stable latitudes. This point 
is not entirely certain, however, and may warrant further anal- 
ysis. 
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The stabilizing effect of the stratification is thus re- 
duced by thermal diffusion, and instabilities take place 
more readily. This has apparently been realized first by 
Townsend (1958), who used the argument in the context 
of shear flows in a stratified atmosphere. It has been ex- 
tended to the astrophysical case by Zahn (1974, 1983). 
In the stellar evolution community, Zahn's expression for 
shear instability in a low-Prandlt number environment is 
sometimes referred to by the curious name 'secular shear 
instability'. 

Applying the argument to the present case of Tayler 
instability, where the stratification also plays a dominant 
role, we repeat the heuristic argument in the preceding 
section by substituting iVc for N in (||). We simplify the 
algebra by assuming tt ^ t. This leads to an erroneously 
high value for iVc in cases where tt happens to be larger 
than the instability time scale r. We correct for this at 
the end of the analysis in Sect. IJ. The minimum radial 
wavenumber kc at which the instability can take place is 
then 



(fcr) 



(10) 



The corresponding length scale 1^ = is the largest ra- 
dial length scale on which the instability can take place. 
Note that it is now independent of the strength of the az- 
imuthal field, contrary to the case when thermal diffusion 
is neglected (eq ||). The instability condition becomes, by 
the same derivation as before: 



(11) 



This is the condition derived in Spruit (1999)^. The con- 
ditions for validity are 



(12) 



and the buoyancy due to composition gradients has to be 
negligible (a different expression holds if the composition 
gradient dominates the buoyancy). Since tj/k <^ 1, condi- 
tion ([ll| ) is less restrictive than (||), and instability sets in 
at lower field strengths. 

Both conditions (||) and (|ll|), are approximate and ig- 
nore multiplying factors of order unity. More exact con- 
ditions can be derived from Acheson's dispersion relation 
(see Appendix in Spruit 1999). 

2.5. Field generation 

2.5.1. The field-amplification loop 

Consider the dynamo process as starting with the Tayler 
instability of an azimuthal field, which itself is produced 
by the winding-up of a seed field. The instability gener- 
ates a new field component whose length scale in the ra- 
dial direction is small because of the strong effect of the 

^ Note that there is a typo in Eq. (49) of that paper. The last 
occurrence of N should be an fi. The correct result is given in 
the appendix of that paper in Eq. (A29). 



stratification (cf eq 10). This instability-generated, small 
scale field has zero average. The differential rotation acts 
on this field, winding it up into a new contribution to 
the azimuthal field. This again is unstable, thus closing 
the dynamo loop. Once the dynamo process has built up, 
the original seed field is unimportant. All three magnetic 
field components then have small length scales in the r- 
direction, but much larger scales in the horizontal direc- 
tions. 

The energy for the dynamo process is fed in by differ- 
ential rotation only, and the small scale field is generated 
by instability of the azimuthal field component itself. This 
is fundamentally different from a dynamo driven by an im- 
posed small scale velocity field such as convection. 

The process of field generation can be broken down 
into two conceptual steps. In the first step, we ask our- 
selves at what level the Tayler instability will saturate, 
if the strength of the initial azimuthal field is given. This 
level determines the amplitude of the small scale processes 
generated by the instability. These in turn determine the 
rate of decay of a given azimuthal field by these processes. 

In the second step, we ask how fast the azimuthal field 
is regenerated by the differential rotation. The process 
here is simply the stretching of the small scale magnetic 
field by the differential rotation. The argument is then 
wrapped up by requiring the regeneration rate to match 
the decay rate. This yields the equilibrium field strength 
of the dynamo process. 

2.5.2. Analogy: convection in a stellar envelope 

These steps can be illustrated with the analogy of stel- 
lar convection. The Rayleigh number is Ra = gH^{'\/ — 
Va)/(j^K), where v and k are the viscosity and thermal 
diffusivity as before, and g the acceleration of gravity. 
Instability sets in when Ra reaches the critical value Rac- 
The typical convective velocity and length scale Vc and 
Ic can be thought of as acting like an effective diffusiv- 
ity Vc ~ vJc- For Ra ^ Rac, we now estimate the am- 
plitude of convection by assuming that the effective dif- 
fusivity Vc ~ f^c is just so large that Ra = Ra^ when 
the effective diffusivities are used in the expression for 
Ra. In other words, the diffusivity becomes so large that 
the effective Rayleigh number is just the critical value for 
onset of instability. With a typical length scale Ic = H, 
one verifies that this yields a convective energy flux of 
Fc = p{gH)^/'^{\/ - Va)^''^. Up to a numerical factor of 
order unity, this is the mixing length expression for the 
convective flux. This well-known 'effective Rayleigh num- 
ber' argument is thus equivalent to a mixing length esti- 
mate. 

While the previous argument yields the convective flux 
when the entropy gradient V — Va is given, this flux is 
usually flxed by other factors. Requiring the flux to equal 
the value that follows from the star's luminosity then de- 
termines the value of V — Va, as well as the convective 
velocity amplitude Vc- 
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We now apply the analogous steps to the magnetic 
dynamo process. In this analogy, the differential rota- 
tion plays the role of the energy flux (assumed to be im- 
posed). The azimuthal magnetic field strength produced 
by winding-up drives the instability and corresponds to 
the entropy gradient, while the other magnetic field com- 
ponents produced by Tayler instability play the role of the 
convective velocity field. 

2.5.3. Instability amplitude for given azimuthal field 
strength 

If the field strength is sufficiently above the critical field 
for Tayler instability there is a range of radial length scales 
on which it can operate. The largest of these is the length 
scale on which the stratification suppresses instability, 1^ 
(eqs The smallest is the scale on which magnetic 

diffusion suppresses instability. Between these scales, the 
growth rate of the instability is not affected much by either 
of these factors and is therefore of the order a = 
With the induction equation, we find that the radial field 
component produced from the initial azimuthal field by 
the unstable displacements is 



Br « B^l/r. 



(13) 



It grows on the instability time scale cr"^, and is largest 
for the maximal radial length scale ^c. As in the convec- 
tive example above we now associate an effective 'turbu- 
lent' diffusivity rje with the small scale field. The value of 
rje is then given, in the same way, by assuming marginal 
stability for the given azimuthal field strength, when rj is 
replaced by rj^. From the stability conditions (||) and ( pi] ) 
this yields 



rjco = r S2( — ) ( — ) 



n 



and 

Vol = 



fl. 



(14) 



(15) 



respectively for the case when thermal diffusion is unim- 
portant (i.e. when the effects of stratification is dominated 
by the composition gradient) and when it dominates. 

One verifies that these expressions are equivalent to 
setting r]e ~ l^cr, with 1^ taken from (10) and (||). As in 
the case of convection, the derivation in terms of marginal 
stability with effective diffusivities is thus equivalent to a 
simple length-and-time scale argument, but perhaps phys- 
ically more illuminating. 

Because of the dominant effects of rotation and the sta- 
ble stratification, all processes are highly anisotropic. The 
effective diffusivity r]e is to be interpreted as the relevant 
diffusivity for gradients in the radial direction. 

2.5.4. Amplitude of the dynamo-generated field 

The field strength in the above is still arbitrary. For high 
field strengths, the effective diffusivity is large, since it is 



driven by magnetic instability, and thus the decay time 
of the magnetic field by the instability-driven small scale 
processes is short. The rate at which the field is amplified, 
on the other hand, is given by the imposed differential 
rotation. For a steady equilibrium, the damping time scale 
has to match the amplification time scale. This condition 
sets the strength of the magnetic field at which the process 
saturates. Since we are assuming that the rotation rate is 
a function of r only ('shellular' rotation), only the radial 
field component contributes to amplification by field line 
stretching. The time scale Ta on which Br is amplified into 
an azimuthal field of the same order as the already existing 
B^ is given by 



Ts.rdrVl ^ B^/Br 



(16) 



The shortest amplification time scale is obtained for the 
largest Br, which occurs for the largest radial length 
scale on which instability takes place. Hence, with ([l3|), 
B^/Br — r/lc- The damping time scale at the same radial 
length scale is 

Td ~ lllVo ~ =. Q/loI. (17) 
Equating the two time scales Ta,Td yields 



when thermal diffusion can be ignored, and 

when it dominates. Here q is the dimensionless differential 
rotation rate 

rdr^ 



(18) 



(19) 



q = 



(20) 



The corresponding azimuthal and radial field strengths are 
B^ — r(47r/9)^/^WA, Br — B^/{kr), hence 



B,,^r{A.pY/^qn^/N, |^=g(§)^ 



(21) 



B,,^r{A^pf/^^qy\'-^f/\-^f'\ (22) 



l/2(-_[^^l/8/ ^l/ 



(23) 



for the two cases, respectively. These are the field strength 
we expect to be generated by the dynamo process. 

2.6. Conditions for the dynamo to operate 

In the derivation I have assumed that the differential rota- 
tion is strong enough to maintain a dynamo process. For 
this to be the case, the azimuthal field produced must be 
large enough for the basic ingredient, Tayler instability, to 
operate. The field produced (eqs |2l| , p^ ) can be compared 
to the minimum field Be for Tayler instability (see |q,pT|): 



Bo_ 

Bco 



^Ac 



(— ) 



1/4 



(24) 
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Bel 



V 



1/8 



(25) 



Setting B/Bc — 1 then yields the minimal shear rate 
needed for the process to operate: 



n - C^^7/4/ V ^l/4 

90 -(^) (^) , 



)l/4(!^)3/4^ 



(26) 



(27) 



( |31| , p2D a re to be interpreted in the sense of their defini- 
tion ~^y. they deter mine the radial angular momentum 
transport under the assumption of 'shellular' (latitude- 
independent) rotation. Since the magnetic stresses are 
much larger in the horizontal directions than in the ra- 
dial, it is likely that they will be quite effective at wiping 
out horizontal gradients in Q, (as suggested in fact by the 
observed rotation in the solar core), so that this assump- 
tion is justified except close to a differentially rotating 
convective envelope. 



2. 7. angular momentum transport 

The main interest of the magnetic fields produced by the 
dynamo process described lies in their ability to trans- 
port angular momentum. Since the azimuthal field in the 
process is built up through the stretching of radial fields 
by differential rotation, the Br and i?^ components are 
nearly maximally correlated. Thus the azimuthal stress 
due to the field generated by the dynamo is, for the two 
cases: 



BroB^O 2 2 3/^\4 

So^^^=pn r g (-) , 



5,«pf7V,(^)V2(-A_) 



1/2 



(29) 



Comparing these expressions, there appear to be cases 
where the stress evaluated with (^9|) is lower than in case 
(p8|), in which the effect of thermal diffusion is neglected. 
This is contrary to expectation, since thermal diffusion fa- 
cilitates the instability process. This is an artefact, result- 
ing from the algebraic simplification made in Sect. 2.5.3, 
where we have assumed that tt ^ t in the expression 
/ = 1 -I- t/tt , so that the effective buoyancy frequency 
is much less than N . When tt > t, the effect of thermal 
diffusion is absent, and instead one should use Nc — N. 
This can be corrected for, approximately, by replacing ( p9[ ) 
by ( psf ) when the latter is larger (see below). 

The magnetic torque on a sphere at radius r is obtained 
from S as usual by multiplying with the lever arm r sin 6 
and integrating over the sphere. The stress can also be 
written in terms of an effective viscosity v^, through the 
relation 



S — pv^rdr^. 
This yields: 



^^^=^2^(il)l/2( « )l/2 



(30) 



(31) 



(32) 



In case 1, the effective viscosity happens to be inde- 
pendent of the shear rate, like a real Newtonian viscosity. 
This is somewhat of a coincidence, as the comparison with 
case shows. 

It should be stressed that the effective viscosity intro- 
duced in this way is extremely anisotropic. Expressions 



3. Applications 

For practical application in stellar interiors, the results 
have to be expressed in terms of the adiabatic buoy- 
ancy frequency TV = g[dp/dP — (dp/dP)ad]^^^ and the 
'compositional' buoyancy frequency = (gdln/i/dz)"'^/^. 
Define the thermal part iVx of the buoyancy frequency by 
7V| = - N^. The analysis of Sect. |^ can then in 
principle be repeated by using as an effective buoyancy 
frequency: 



(28) N^^^N'^+N^/il + r/rr), 



(33) 



where r is again the instability time scale and tt the ther- 
mal diffusion time scale on length scale I. This takes into 
account that the thermal part of the restoring force is re- 
laxed on long time scales and short length scales, while 
the compositional part can not be reduced. 

3.1. Effective viscosity 

This can be carried through to obtain the generalization 
of expressions (^ 29) for the stress. The algebraic com- 
plexity of the resulting expressions, however, is incom- 
mensurate with the level of sophistication of the present 
analysis. The dependences on N^^ and are probably 
monotonic, so that it is simpler to use a patching for- 
mula to connect the two limiting cases. Taking also into 
account that the dynamo process requires a minimum ro- 
tation gradient qaiin to operate, a possible patching for- 
mula for the effective radial viscosity Vre produced by the 
dynamo-generated magnetic field is: 



where 



r.ei = r2fimax[(-^)i/2( 



and 







[q < gmin)- 



(34) 



(35) 



(36) 



(37) 



The factor / causes the stress to vanish smoothly as the 
gradient of the rotation rate approaches the minimum 
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value required for dynamo action. The max operator in 
( ^ ) corrects for the error made in deriving ( 29 ) for cases 
where thermal diffusion has no effect, as discussed in 2.7. 
The first factor in 



( |34| ) guarantees that the stabilizing ef- 
fects of both the compositional and the thermal stratifi- 
cation on the instability process are taken into account in 
the final result. 

For gmin one can use a patching formula to connect eqs 
(p^,p^, for example 



90 + 91, 



where 



,N„ 



^0 = 



91 - {^?'\-P^f''C-f''. 



(38) 

(39) 
(40) 



The sum in ( pq ) takes into account that dynamo action is 
possible only when the rotation gradient is strong enough 
to overcome both the stabilizing effects of the /i— gradient 
and of the thermal buoyancy. 

3.2. Mixing 

The fluid motions involved in the dynamo process also 
imply a certain amount of mixing. The mixing in the ra- 
dial direction is effected by the same displacements that 
produce the effective magnetic diffusivity. Up to a factor 
of unity we can set the effective diffusivity De for mixing 
in gradients of composition equal to rjc. With (1^15|) and 
( |l8| , pj| ), patching in the same way as Eqs. (^^) yields 
for the effective diffusivity, 



■fill 



where / is defined in 

n 



and 



Del = r^r!max[g(-^)3/4(__)^/^g4( " ) 



^3/4 



n 



(41) 



(42) 



(43) 



One verifies that the max operator in ( [43|) switches be- 
tween its arguments under the same conditions as it does 
in (H). 

From these expressions one finds that Dc is much 
smaller than v^, approaching it only when r/Ndil./dr is 
of order unity. This reflects the fact that the angular mo- 
mentum transport is done by magnetic stresses, which are 
much more effective than the Reynolds stresses associated 
with the fluid motions in the dynamo process. 

4. Discussion of the result 

The magnetic stress or its equivalent effective viscosity 
( ^ H^) is the net effect of several processes: the winding- 
up of radial fields by differential rotation, the instability 



of the resulting azimuthal field, and the internal dissipa- 
tion of small scale fields by reconnection. Together these 
constitute a dynamo process which generates a magnetic 
field whose strength is governed by the differential rota- 
tion. Since the length scale of the field turns out to be 
small in the radial direction, the radial transport of angu- 
lar momentum by the magnetic stresses can be represented 
by an effective radial viscosity [as defined in (34)]. Both 
the length scales and the fleld strength itself are much 
larger in the horizontal directions {9, (f>). Variations in the 
rotation rate on horizontal surfaces are therefore likely to 
be quite small. 

Comparison of ( p5| ) and (^) shows that the presence 
of a gradient in composition (/i) has a dramatic effect. As 
long as this gradient can be neglected the effective radial 
viscosity scales as £7"^/^ (for fixed dimensionless differential 
rotation q); when it dominates, the dependence is as fl^. 
The steep dependence on can be accounted for by noting 
that the stress S is proportional to Br and B^, while both 
B^ and Br/B^ contain a factor {fl/N)'^. The stabilizing 
effect of the stratification thus effectively appears three 
times. 

The dependence on Q is steeper than in some purely 
hydrodynamic transport mechanisms. On the other hand, 
the scale of the effective viscosity multiplying this depen- 
dence, of the order r^flq^, is much larger. To see what 
this implies, consider for comparison Zahn's (1973, 1983) 
formalism for shear instability in the presence of thermal 
diffusion. Assuming a critical Reynolds number Rcc ~ 10^ 
for small scale turbulence, this process yields an effective 
viscosity of 



(44) 



(45) 



Comparing gmin.z with (|27l), one finds that the dynamo 
process has the lower value for the critical shear rate when 



while the critical shear rate for it to work is 

9min,Z ^ -^{RSc-) ' ■ 



i-r 



> 



V 



(46) 



Since v/rj ^ 1, and K/{r'^N) is a very small number, the 
dynamo process generates an effective viscosity already at 
much lower rotation rates. In the present Sun, for exam- 
ple it operates, while hydrodynamic turbulence according 
to Zahn's formalism is excluded by condition (p5[). Under 
conditions when both hydrodynamic turbulence according 
to Zahn's estimate and the dynamo process operate, the 
effective viscosity of the dynamo process is larger, pro- 
vided that [compare (H,!!)]: 



- ( 



N 



1/3 



(47) 



For solar parameters, (^/(r^A^))^/^ « 10^'^, hence the ef- 
fective viscosity of the dynamo process is also larger, when 
both processes are possible. 
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The high effectiveness of the dynamo process in the ab- 
sence of composition gradients is likely to remove most of 
the differential rotation in compositionally homogeneous 
layers of a star, with all the gradients concentrating into 
the inhomogeneous layers. The steep dependence on f2 will 
cause a certain 'convergence': somewhat independent of 
other factors influencing the rotation, there will be a ten- 
dency for il/Nf^ to evolve into a limited range. 

The process as described here makes the essential as- 
sumption that the initial seed fleld being wound up is suf- 
ficiently weak. For strong initial fields, such as those ob- 
served in magnetic A-stars and white dwarfs, the process 
does not apply. Instead, in such stars rotation is likely to 
become and remain homogenous on a short time scale (cf 
discussion in Spruit 1999). 

Assuming that the initial field is indeed weak, it takes 
a while for the dynamo process to reach its equilibrium 
amplitude. This may still be a sufficiently short time scale 
in most circumstances. 

Since the dynamo process produces its own 'turbu- 
lent diffusion', it enhances the decay of an initial field|^. 
Thus the process also enhances the distinction between 
two regimes: an initially weak field serves as a seed field 
but is eliminated by the action of the dynamo process, 
while a stronger initial field will prevent the dynamo pro- 
cess from operating altogether. 

5. Estimates for the Sun 

For the present Sun, with a rotation rate = 3 • 10~^, 
the buoyancy frequency in the outer core, with radius 
r = 5 • IQi", is of the order N = IQ-^, while /t « 2 • 10^, 
77 w 1.5 • 10'^. In the so-called tachocline (Spiegel & Zahn 
1992), the rotation rate changes from a pattern with a 
variation of 30% between equator and pole in the convec- 
tion zone, to a uniform rotation in the core. The width 
of this transition is believed (Charbonneau et al. 1999) to 
be of the order O.OSi?©, so that the shear rate is g ~ 2 
in the tachocline (but varying and changing sign with lat- 
itude). In this region of the core, composition gradients 
(mostly due to the gravitational settling of Helium) are 
quite small, so that case 1 apphes. With these numbers, 
we find 

^ « 50. (48) 

so that the conditions for dynamo action are satisfied. The 
expected field strengths are 

« 1.5 ■ W G, Br^l G. (49) 

With this field strength, the Alfven frequency is uja ~ 
10^^ s^^. The actual frequency of Alfvenic modes is re- 
duced by the strong Coriolis force. Waves with length scale 
~ r have a frequency 

cr = Lul/n ~ 3 • 10"^ (50) 

Strictly speaking: only the non-potential parts of this field. 
If the dynamo operates only in a limited range of radii, it can 
not eliminate a field that is anchored in other parts of the star 



which is also the characteristic growth rate of the insta- 
bility. This corresponds to a time scale of the order 100 
years. 

6. Discussion 

As in the case of accretion disks, an externally imposed 
small scale velocity field is not necessary for obtaining field 
amplification by a dynamo process in stars. In addition to 
differential rotation, the only ingredient required is an in- 
ternal instability in a toroidal magnetic field, to provide 
new poloidal field components on which differential rota- 
tion can act. Several distinct instabilities are candidates 
for this key ingredient: Balbus-Hawley instability, buoy- 
ant (Parker-) instability, and the toroidal-field instability 
studied by Tayler. I have argued here and in Spruit (1999) 
that the latter is probably the most relevant, and have 
shown how a small-scale magnetic field may be produced 
in stably stratified layers in stars through this instability, 
with differential rotation as the energy source of the pro- 
cess. The required differential rotation may be produced 
by an external torque like a stellar wind, or by the internal 
evolution of the star. The main significance of the process 
thus lies in its effect on the internal rotation of the star, 
and for questions such as the amount of rotation to be 
expected in the cores of giant stars. This in turn is critical 
for the interpretation of the rotation of pulsars (Spruit & 
Phinney 1998) and white dwarfs (Spruit 1998). 

The horizontal components of the magnetic field pro- 
duced by this process are much larger than the radial com- 
ponent, and the radial length scale of the field is much 
smaller than the horizontal length scales. Because of the 
small radial length scale, the effects of the microscopic 
(atomic) magnetic and thermal diffusion coefficients have 
to be included in the analysis. Because of these effects the 
gradient in the rotation rate must exceed a certain mini- 
mum for dynamo action to set in. The measured gradients 
in the tachocline of the present, slowly rotating Sun are 
large enough, but only by a factor of 100 or so. Because of 
the large ratio of horizontal to radial field strengths, the 
angular momentum transport by the magnetic stresses is 
very anisotropic, and far more effective at smoothing out 
horizontal variations in rotation than radial variations. 
This is interpreted here as the main reason for the pe- 
culiar abrupt transition from the rotation pattern in the 
convection zone to that in the core of the Sun. 

It must be stressed that the quantitative estimates of 
the dynamo process made here ignore all multiplying fac- 
tors of order unity. Experience in similar cases shows that 
these factors can sometimes compound to rather large 
numbers. Also, I have essentially ignored the difference be- 
tween the spherical and the cylindrical radial coordinate 
in making estimates. At the level of the present analysis, 
there is no justification for greater quantitative detail, but 
the resulting uncertainty must be borne in mind. Finally, 
the analysis applies only for cases where the initial mag- 
netic field of the star is small enough. In strong magnetic 



10 



H.C. Spruit: Dynamo action in a stellar interior 



fields like those of the magnetic A-stars and white dwarfs, 
the dynamo process is suppressed. 
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